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Abstract 

We derive a low-energy Hamiltonian for the elastic energy of a Neel domain 
wall in a thin film with in-plane magnetization, where we consider the con- 
tribution of the long-range dipolar interaction beyond the quadratic approxi- 
mation. We show that such a Hamiltonian is analogous to the Hamiltonian of 
a one-dimensional polaron in an external random potential. We use a replica 
variational method to compute the roughening exponent of the domain wall 

for the case of two-dimensional dipolar interactions. 
05.70.Np, 68.55.Ln 
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Typeset using REVT^ 



I. INTRODUCTION 



The roughening properties of elastic manifolds in the presence of quenched disorder is a 
well studied problem. The subject of study is usually the asymptotic (long-distance) regime 
in which one wishes to compute the roughening exponent ( of the D-dimensional elastic 
manifold (described by a displacement field (p{x)) defined by 



mx)~<f>{x')y)=\x-x'\'<, (1) 

where x is the position vector in the D-dimensional configuration space and where the 
brackets represent an average over thermal fluctuations and the overline an average over 
the realizations of the quenched disorder. Obviously, the study of the crossover between 
short distances, where the thermal fluctuations dominate and the long-distance behaviour, 
which is determined by the disorder, is also of interest. Among the methods applied to the 
study of this problem are the use of Imry-Ma type of argumentsEHi, the mapping of one- 
dimensional interfaces in two dimensions to the noisy Burgers equation^, which yields an 
exact result for the roughening exponentii, variational approaches involving replica averag- 
ing and replica symmetry breakingi^ and functional renormalization group calculations0'0. 
The role played by long-range interactions does not seem to be so well understood, although 
it follows from the one-loop functional RG results of Emig and NattermannS that in the 
important case of a magnetic domain-wall in the presence of long-range dipolar interactions, 
one simply needs to replace the expansion parameter e = 4 — D hj e = 3{3 — D) /2 (where D 
is the dimension of the manifold) in the expressions for the critical exponents obtained in the 
absence of dipolar interactions, in order to account for the presence of these interact ionsEl. 
In particular, if one applies this result to the case of a line domain in a thin magnetic film 



{D = 1), the effective expansion parameter is e = 3 in both cases, which means that in 
this case and to one-loop order, the critical exponents obtained are the same in the absence 
or presence of dipolar interactions. Although one expects the RG calculations to provide 
a qualitative understanding at low dimensions compared to the upper critical dimension 
= 3 in this case), it is still questionable to what extend one can perform such ana- 
lytical continuation using just the one-loop results. Furthermore, the results of Emig and 
Nattermann were based on a series expansion of the dipolar energ which meant that, 
in the absence of disorder, the Hamiltonian describing the low energy degrees of freedom 
of the domain wall was quadratic in the domain-wall displacement field, i.e. the random 
potential is the only source of non-linearity in the problem. One still needs to justify that 
such approximation is sufficientlll. 

Having in mind such difficulties, we wish to discuss the properties of a one-dimensional 
domain wall in a two-dimensional ferromagnetic film, in the presence of dipolar interactions 
and a short-range correlated random field, where the dipolar interaction is treated beyond the 
quadratic approximation. We wish to check whether it is still possible to map such a problem 
to a (modified) Burgers equation, as was done by Huse et al.li for the case of random-bond 
disorder and by Zhangi for the case of random-field disorder, in the absence of long-range 
interactions. Such a mapping was achieved because one can map the line domain to the 
world line of a quantum particle in 1+1 dimensions. The free energy of the line domain 
can then be shown to obey a Schrodinger equation in imaginary time in the presence of a 
random potential. A Cole-Hopf transformation can then be used to convert this equation 
into a noisy Burgers equation with additive noise, for which the stationary distribution is 
known in one-dimensionS. This allows one to compute the roughening exponent. In 



our 



case, the presence of the dipolar interaction does not allow one to write the energy of the 
domain wall as an action for a single quantum particle, but instead we obtain an action for 
a quantum particle which interacts with a annealed 'vector potential', whose propagator is 
determined by the dipolar propagator. In 1+1 dimensions, a gauge transformation allows 
us to substitute the vector potential by a scalar one, which can then be integrated out, 
leaving us with an action similar to the action of a polaron as considered by FeynmanllBlli, 
i.e. the dipolar interaction introduces a 'self-retarded' interaction of the particle with itself. 
Despite the fact that such action cannot be integrated exactly, we can still apply a variational 
method in the spirit of FeynmanS^, in which one introduces replicas of the system in order 
to average the free energy of the domain wall over the realizations of the random field. Using 
the hierarchical replica symmetry breaking ansatz of Parisi0, we were able to derive a set 
of non-linear equations whose self-consistent solution yields the long-distance physics of the 
problem within the realm of this approximation. With this method, we are able to treat 
two different situations. The one of a thin ferromagnetic film grown on a non-magnetic 
substrate, in which the dipolar interaction has the usual three dimensional form, and the 
one of a thin ferromagnetic film grown between two type I superconductors, i.e. the case 
of a superconductor /ferromagnet /superconductor (SC/FM/SC) heterostructure. We have 
shown elsewhere^ that in this case the dipolar interaction has a two-dimensional form at 
small wave-vectors, coming from the renormalization of the magnetic energy due to the 
Meissner effect in the superconductors. Although we were unable to obtain an analytic 
solution of the variational non-linear equations, we were able to justify that in the case of 
two-dimensional dipolar interactions (i.e. for the SC/FM/SC heterostructure), the quadratic 
approximation is indeed adequate. This case is not only easier to treat analytically, it is 



also more interesting than the case of three dimensional dipolar interactions, for which 
one expects the dipolar interaction to be irrelevant (see above). Therefore, we restrict our 
discussion to the case of two-dimensional dipolar interactions in this paper. Furthermore, we 
can also show that in this case, if one works within the quadratic approximation, the elastic 
energy of the domain- wall has the same analytical form as the energy of a domain wall in a 
ferroelastic material as studied by Kolomeiskii et al.0 and also of the energy of a liquid-gas 
interface in a disordered solid, as studied by Hazareesing and Mezardi. In both cases, one 
obtains, in two dimensions, a value of 1/3 for the roughening exponent. This is indeed the 
value that we obtain for the roughening exponent in our problem. It should be nevertheless 
noticed that these three problems, though mathematically similar have a different physical 
origin and motivation. 

The structure of this paper is as follows: in section II, we derive a low-energy Hamil- 
tonian describing the Neel domain-wall elastic degrees of freedom. This Hamiltonian can 
be interpreted as the action in imaginary time of a quantum particle in 1+1 dimensions. 
In section III, we use the replica trick to average over the random field realizations, which 
generates a many particle non-quadratic action. In section IV, we introduce a generaliza- 
tion of the variational ansatz of Mezard and Parisi which allows us to consider the effect of 
the short-range correlated disorder and of the non-local dipolar term in the many particle 
action. One obtains a set of non-linear equations which have to be solved self-consistently. 
In section V, we perform the analysis of these equations in the particular case of the hi- 
erarchical replica symmetry breaking solution of Parisi and we discuss the validity of the 
quadratic approximation for the dipolar term in these equations, which permits to deter- 
mine the roughening exponent of the domain wall. Finally, in section VI, we present our 
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conclusions. 



II. ENERGY OF A SINGLE DISTORTED DOMAIN WALL IN A 
TWO-DIMENSIONAL XY MODEL WITH DIPOLAR INTERACTIONS AND IN 

A RANDOM FIELD 

We consider a Neel (i.e. 180 degrees) domain-wall in a thin two-dimensional ferromag- 
netic film in the presence of dipolar interactions and a short-range correlated quenched 
random field. The magnetization is in the plane of the film. In the straight domain confi- 
guration, the wall is oriented along the y axis and the center of the wall is located at x = 0. 
The thermal and random field fluctuations give rise to deviations of the domain-wall from 
its straight configuration, i.e. at a given position y, the center of the domain is displaced 
from X = to (we neglect overhangs). In this case, equation (|I]) reduces to 

{{x{y)-x{y')Y)=\y-y'\'^ . (2) 

The distance \ y — y' \is considered to be very large with respect to the lattice spacing and 
to the length over which the random field is correlated (see below). In the absence of dipolar 
interactions, one can, as stated above, through a mapping to a Schrodinger equation§, show 
that C = 1- 

In order to determine the elastic energy of the domain-wall, we consider the underlying 
system of ferromagnetic spins to be in a square lattice with lattice constant a and overall 
dimensions LxX Ly. The system has anti-periodic boundary conditions along the x direction 
(in order that there is a single domain-wall in the system) and periodic boundary conditions 
along the y direction and is at finite temperature T. Such system is described by a (classical) 
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XY model Hamiltonian, with dipolar interactions and in a random field, 

= — - ^ Jij Si - Sj — fi ^ hi ■ Si (3) 

<i,j> i 

\^ f ' _ ^ {Rj - Rj) ■ Si {Rj - Rj) ■ Sj \ _ T-) / C-S/N2 

2 ^\\R,-R,\^ \R,- R. h+2 J 2^^ ^> ■ 

i-^j \ I J I I J I / J 

For the case of two-dimensional dipolar interactions (i.e. for the SC/FM/SC heterostruc- 
ture), the constants and G are given in the MKS system by /i = iii^\iB9L and G = 
ath^Xl ^dLfJ'BY, where /io is the magnetic permitivity of the vacuum, fiB the Bohr magneton 
and ql the Lande factor of the spin system. The constants and fir are, respectively, 
the London penetration depth and relative magnetic permitivity of the superconductor. In 
this case, n = 2. In the case of three dimensional dipolar interactions, which we will not 
consider, G = ^{fis gi)"^ and n is equal to 3. The quantities Jij and D have the dimensions 
of energyi^. The field hi is a random field with Gaussian correlations, i.e. h" hj = Aij6"^. 
In order to stabilize the domain-wall, we have added to the Hamiltonian above an easy-axis 
anisotropy term, along the y axis, represented by the last term in (^. 
For a straight Neel domain-wall, the spin configuration has the form 

= S cos <f){xi) (4) 
Sf = S sin (5) 

where i is the site label and where (/>(x) goes from —7t/2 when x — > —L^/2 to tt/2 when 
X Lx/2, with S being the value of the spin. The function changes its value from 
— 7r/2to7r/2ina region of width w around x = 0. Since w is typically of the order of tens of 
nanometerJ^i and we are only interested in the long-distance properties of the domain-wall, 
we can write = 7r5{x). If the center of the domain- wall is displaced from = to 
Xi = x{yi), the spin configuration becomes 



Si = S COS (j){xi - x{yi)) (6) 
Sf = S sin (j){xi - x{yi)) . (7) 

We wish to compute the partition function Z = Tr{e~^'^^''^'^), where SH = SHexc + SHrf + 
STi-dip is the energy difference between a configuration described by the displacement x{y) 
and the straight domain configuration. We define the discrete Fourier transforms 

Ri 

where q belongs to the first Brillouin zone. We will approximate the sums over i by integrals, 
according to the formulas, 

valid for a square lattice. It can be shown that the (positive) contribution to the exchange 
energy coming from the domain-wall deformation is given to second order in the displacement 
field, byi 

^Hexc = -^1^2J('^(^^'^i') ~ Ji(lx,0))x{qy)x{-qy) 



q 



2Lya 

1y 



^qlx{qy)x{-qy) (11) 



where N = LxLy/a"^ is the total number of sites and where we have expanded J{q) ~ 
Jo — Ja?q^- In real space, equation (|TT|) reads 

which has the same form as the kinetic part of the action of a particle of mass a in one 
dimension, with a = Ti'^JS^/a being the line-tension of the domain wall. The contribution 



coming from the random field can be sliown to be equal tc 



2 



STirf = I dy I dxh^{x,y) (13) 



2 



where fl = fJ,B giT^S/ a and = / a0. This formula can be made plausible if we consider 
instead a constant field along y. It can also be easily seen from (|TI]) that the easy-axis 
anisotropy term does not contribute to the deformation energy of the domain-wall. 

Since we wish to compute the partition function, we have to consider the factor 
g-Hdip/fesT^ where T-tdip is given by 

^ q 

where V = L^Ly is the area of the sample. For two-dimensional dipolar interactions, g = 
2^f^(fi'L/^B)^ and J-'(q) = For three dimensional dipolar interactions, g = ^{gL^^BY 

and J^{q) = ^. Thus, by considering the different dipolar kernels J-'{q), one can treat the 
two different cases. 

Using a Hubbard-Stratonovich transformation, we write e~^''»p/'^^-^ as 



(15) 



where is a normalization factor. We can now show that 



Ly 



q.Sq = '-^ r dye-''^^^^y^-"^yy^, (16) 
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from which it follows that 



q "-- 
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and one can see that there is no contribution to the dipolar energy when x{y) = const. The 
partition function of the domain wall's degrees of freedom involves the functional integration 
of the Boltzmann factor e"''^/*^-^"^, over the space of functions x{y). 
We shall now make a change in our notation. We will write 

y^T, x{y) x{t), ^ x(r), 

dy 

r^SJg 

Qy ^ujm, e = ■ (18) 

The subscript m stands for the fact that the wave- vectors uim are discrete, due to the periodic 
boundary conditions along r {y). If we use this notation in the formulae above, we have, 
collecting all the factors, the following expression for the partition function 

Z=M j Vxir)VA{q,,Um)e-^/^''^ (19) 

where the action S is given by 

S = -o I dTx{T)+fi dr dxh^{x,T) 

2 J-^q- J-if Jo 

r¥ 1 

-ie dr A{x{r),T) x{r) + — A{q^,LJrn) ^'^{qx,^^m) A{-q^, -Urn) ■ (20) 
J_hi zV 

This is the action of a single quantum particle in 1+1 dimensions, with 'mass' a and 'charge' 
e, at 'inverse temperature' Ly and with a 'Planck constant' equal to k^T . Such particle is in a 
external random potential V(x, t) = [i dx' h^^x', r) and interacts with a 'vector potential' 
)i characterized by a propagator equal to J-'{qx,LJm)- This result belongs to the 
well known class of mappings of 2 dimensional problems of classical statistical mechanics 
to 1 + 1 dimensional quantum problems. Now, since one has periodic boundary conditions 
along the r {y) axis, the following identity holds 

dr A{x{t), r) x{t) = - [ ' dr T^"^ dx' ^"^^^'^ , (21) 
y. .I^ht Jo OT 
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which corresponds to a gauge transformation of the action!!. The use of this identity ehm- 
inates the dependence on the 'velocity' x(r) in all terms of equation (^) except the first 
one. The action S now reads 

S = \al' dTx\T)+ dTV{x{T),T)+ie [ ' dT r " dx' ^^^^ ' 



tji I — kn. In 

2 2 " 



Or 

+ ^ 5Z Mlx^^m) ^~^{qx,^^m) A{-q^,-Um) ■ (22) 

We can now integrate out the field A^q^^Urn)-, obtaining an action entirely in terms of x(r) 
alone. We obtain 

Ly Ly 

S = -a[' dTx'^{r)+ [' rfr1/(x(r),r) 

ZV at. I —ksL I— hi 

This action can be used to compute the roughening exponent of the domain wall. Equation 
(p3|) also shows that, apart from the interaction with the r dependent random potential 
V{x, r), the action S is analogous to the action of a one-dimensional polaron, as written by 
FeynmanSlli. 

However, in order to compute the roughening exponent, one has to determine the fol- 
lowing quantity 



/ Vx{t) e-^AsT ^x{r)-x{r')y 



{{x{r) - x{r')f) = ^ Jp^(,);^5a;t (24) 

which involves determining the average of the distribution of quenched disorder over a quo- 
tient. This type of averaging places the same problem as when one needs to average over the 



free energy of the system, i.e. one needs to determine the average \ogZ, instead of simply 
averaging over the partition function Z. This can be seen by including a source field coupled 
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to x{t) in Z. The two-point correlation function appearing in (^) can be computed from 
log 2 by differentiating it with respect to the sources and taking these sources to be zero. 



III. THE AVERAGE OVER DISORDER. REPLICA TRICK. 

We need to perform the average over the random field h^iqx, uJm), which obeys a Gaussian 
distribution with a second moment equal to 

In order to perform the averaging over disorder in the free energy, one introduces n distinct 
copies (replicas) of the system^ and makes use of the identity 

— 1 

logZ = lim . (26) 

n^O n 

One performs the disorder average over for integer n and makes the analytical continu- 
ation to n — > 0. One obtains 



Z"= / JJPx„(r)e-^«/'=«^, (27) 

a 

where the 'replicated' action Spt is given by 

Sn = l^Il n dr{xl{r) + nlxlir)) + -^Y.^^ T cos(g.x.(r)) 

a 2 a,qx 2 

2 2 

2V ^ qi. J_hL J-bL 



and where the last term represents the inter-replica interaction, due to the presence of the 
random potential. We have included in the effect of an applied field hl^^i^x) = (crQl/fi) x. 
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This term gives rise to the second term in the action, which is a harmonic potential. The 
presence of this term is necessary to guarantee that Sr is bounded from below, given that 
the last term of p^) comes with a negative sign. Qq can be set to at the end of the 
calculations. 

IV. A VARIATIONAL ANSATZ FOR THE FREE ENERGY. VARIATIONAL 

EQUATIONS. 

The action (pSD cannot be integrated exactly. Therefore, we need to develop an approx- 
imation scheme. We shall follow Mezard and ParisiB and Goldschmidti. We choose a trial 
quadratic action of the form 

Ly 

So[Va{r)] = laJ2 fl dr{xl{r) + Qlxl{T)) 



a,b 



-E ld^Va{r)x,{r) (29) 

where we have included a source field rjalr), which is useful in the calculation of correlation 
functions. The Mezard-Parisi ansatz has been generalized by including a matrix function 
QabiT—T')-, with a non-trivial dependence on t—t', necessary to take into account the effect of 
the non-local dipolar interaction. In the Mezard and Parisi case, QabiT — t') = Qab S{t — t'), 
since they considered only short-range correlated disorder. Goldschmidt has considered 
several different forms for the kernel Qabi'T — t'), all of which non-local in r — r', and which 
are able to account for the problem of a quantum particle in a disordered static potential 
and for the related problem of a polymer in a random potential with long-range correlated 
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disorder (see below). Since Qabij) = Qabij + Ly), we can represent Qahij) in terms of 
Matsubara modes Qabij) = ^l^y Xlm QS^*'^'"''- The values of fi^ and are to be chosen 
appropriately. We also demand Sq to be translationally invariant when $7^ ^ 0. This gives 
rise to the constraint 

Qlb = -Y.^lb- (30) 

It can be shown that the normalized generating function for the correlation functions of such 
a quadratic action, is given byi^ 

^^[r?a(r)] = exp ['l/2(fcBT)2 J] I f IJt' r]a{r) Gab{r - t') r]t{T')] (31) 

where Gabij — t') = (xa(r)xh(r')) is given in terms of its Matsubara modes, by 

G2 = kBT[a{u:li + Cl') - Q"^]-J, (32) 

where / is the unit matrix in replica space and Cl'^ is the diagonal matrix with elements 
equal to Q^. The free energy as computed from the trial action S'o[0], with the sources rjalr) 
equal to zero, can be shown to be equal to 

Fo = - ^ 5^ In [ {l/ksT){au;l + (aCl' - Q^,) ] (33) 

where /x indicates the basis in which the matrix aVt^ — is diagonal. One now uses the 
Feynman inequality for F = —/c^T log Z", which states that, 

F<F,ar = Fo+{SR-So[0])o, (34) 
where the average in the rhs of (Bl) is over the trial quadratic action 5*0 [0]. Note that 



J-' = —ksT log Z = lim„_+o F/n. One obtains for F^ar the expression 
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cr 

+ 2Z 



2L„ 



(35) 



m,ii 



^ a,?n 



+ 



ksTL 



where the following identities, valid for a quadratic action at zero sources, (e *'?^^'»(^))q = 

^-\qlGaa{0) ^-^^ (^Q-i<lx{Xa{T)-Xi,{Q))\^^ _ g- 1 (Gaa {0) + G66 (0)- 2Ga6 (t)) ^gj^g USed. FollOWlUg 

Hazareesing and Mezardi, we choose the random field to have anisotropic correlations, 
A{(lx-,^m) = A e^^^^'^^'^^, with correlation length B along x and we choose the two- 
dimensional form of the dipolar kernel, i.e. J^{qx,ujm) = 1/(5'^ + ^m)- have to minimize 
(^) with respect to the free parameters and Such procedure gives the following 
self-consistent equations 



q2 _ q2 _ 

"a — "0 



QZ = Qlb + fTSat m ^ 



Ql 



ab 



An 
/27rkBT 



{B' + Gaam 



-1/2 



with 
jl'A 



B' + Gaa{0) + G,,{0)-2Gab{0)) 



-1/2 



^27rkBT 

fa = -e' / ' Ar{ 1 - cos(^„r) 

^0 



r 



e*T""»(^) Erfc 



r 



(36) 
(37) 

(38) 



(39) 
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with Q^^ given by (|30|) and where 7aa(T) = GaaiO) — Gaai^ We have approximated the 
sum over by an integral, a procedure which becomes exact in the thermodynamic hmit, 
and Erfc(x) is the complementary error function. These equations constitute a closed non- 
linear system which has to be solved self-consistently in the limit n — > 0. We consider the 
solution of these equations in this limit in the next section. 

V. ANALYSIS OF THE VARIATIONAL EQUATIONS FOR A 2D DIPOLAR 
INTERACTION. HIERARCHICAL RSB. 

One needs to take the limit n — in the self-consistent equations derived above. This 
implies dealing with the matrix structure of the equations in this limit. One cannot do this 
in general. However, there is a particular parametrization of these matrices, due to ParisiEl, 
which is sufficiently general for our purposes. This is the so called hierarchical RSB. The 
diagonal elements of the matrix Qab (with dimensions n xn where n is arbitrarily large) are 
taken to be all equal to q and the off-diagonal elements are taken to be equal to go- Qab 
is then partioned in block- diagonal submatrices, the elements of Qab outside these blocks 
keeping the value go and the elements inside the blocks taking a new value gi (except the 
diagonal elements which keep their value g) . The procedure is repeated in an equal form for 
every diagonal block submatrix, the off-diagonal elements of these submatrices keeping the 
value gi and the elements inside the diagonal sub-blocks taking the value g2. This procedure 
is repeated 'ad infinitum'. The multiplication rules obeyed by these hierarchical matrices 
are analytically continued to — ^ 0. The off-diagonal elements of Qab are then parametrized 
by a function q{u), u G [0, 1] and the diagonal elements are parametrized by the number g. 
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The constraint 



takes the form 



g = / du q{u) 
'o 



(40) 



and one can show that the matrices are also parametrized by a function g"^{u) and by 
g"^ which parametrizes the diagonal elements. In the thermodynamic limit = Ly ^ oo, 
LOrn —>■ uj, g"^ ^ g{uj) and g'^{u) — > g{uj,u) and one gets, using equations (|36| - ^ ), the 
following relations, where g^r) {g{T,u)) is the Fourier transform of giyj) {g{uj,u)), 



n' = nl--^=^iB^ + ~g{r = o)r/^ 

q{u) = ^^J B' + 2Cg{T = 0) - gir = 0, u)) 



2-nkBT 

POO 

f{uj) = — / (ir ( 1 — cos(u;r) ) 
Jo 



r 



(41) 
(42) 

(43) 



_27r7(r) A^-f^^ir)'^ ' ^^^^ V V/^(^) 
and where 7(r) = g{T = 0) — ^(t). One sees from (|4l| ) that when n ^ 0, Cl'^ = Qq. The 
functions g{uj) and g{u,u) are in turn given by 



q{0) 



^ ,^dv [q]{v) 



g{uj,u) 



v^a{u^ + nl)-f{u) + [q]{v) 

keT ( g(0) r dv 



(44) 



(45) 



«[a(a;2 + fi2)_/(^) + [g](^)] 

where \q\{u) = uq{u) — dvq{v). Substituting (|4^ and (|45|) in equation ( ^2]) and differen- 
tiating with respect to m, one obtains the equations 



q{u) = 



or 



q{u) 



2'K{kBTf'^ J 



2/3 



duj 



(46) 



21X [a{uj^ + nl)-f{u:) + [q]{u)y 



-1/3 



If the first equation is valid in the whole interval u G [0, 1], its solution is simply q{u) = const, 
which is the replica symmetric solution. The second equation is easily solvable if f{uj) = 
(i.e. e = 0) and the solution wa. di.cu.=ed by M&axd and Pari J. Also, one has 

{{x{t) - x{r')y) = hm - J2(Gaa{0) - G^r - r')) 

a 

= 2 7(r-r') (47) 

which allows for the computation of C, once 7(r) is know. In the presence of dipolar in- 
teractions, equation (|4^) depends on 7(r) and one has to resort to a numerical approach 
if one wishes to solve it self-consistently. However, if one simply takes the zeroth-order 
result f{uj) ~ — Y | u |, coming from the expansion of /(cu) in powers of 7(r), one sees 
that, apart from the contribution coming from the line tension a, the equations P^ ) to 
(HSI) are identical to the ones obtained by Hazareesing and Mezardi in their study of the 
roughening properties of a liquid-gas interface in a disordered solid. Also, if one expands 
the last term of equation ( P^D to quadratic order in x{t) — x(r'), one can see more directly 
that the domain-wall energy has, in this order of approximation, the same analytical form 
as that obtained by Kolomeiskii et al. for a domain-wall in a ferroelastic material and as 
that used by Hazareesing and Mezard in the study of the problem referred above (provided 
one neglects the contribution of the line tension in our case). Therefore, if one takes this 
approximation, one can, following Robbins and Joanny and also Kolomeiskii et al. use 
a Imry-Ma type of argumenti to show that the expected value of C is 1/3, this behaviour of 
7(r) being valid for lengths L much larger than the Larkin lengthH ^, which is given in our 
problem by ^ ~ When L ^ ^, we can disregard the line tension contribution to the 

action, since this energy term scales with L^^ . However, at short length scales compared to 
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^, this term cannot be disregarded, since it gives an uJ^ contribution to the denominator of 
the Green's function appearing in equation (^^. This contribution gives rise to a behaviour 
7(r) ~ r, at length scales < i3, as can be seen from the Mezard-Parisi solution@@. Having 
in mind such behaviour of 7(t), we can write (E3|) as 



j{iS) = — e / (ir ( 1 — cos(u;r) ) 
'o 

f oo 

' (ir ( 1 — cosftur 



1 



27r7(r) A:^i^l'^{T\ 
1 r 



T 



27V2(r) 
r 



• (48) 



The integrand in the first term is regular at small r (due to the behaviour of 7(r) at small 
r) and its contribution is unimportant if cu^ ^ 1, which is the set of wave- vectors which 
determines the roughening exponent (. In the second term, we can use, if i3 <^ ^, the 
asymptotic expansion for Erfc(x). We obtain, omiting the regular contribution of the first 



term, 



-,2 /-oo 



dr (1 — cos(ujt) 



67(r) 



+ 



(49) 



where the dots stand for terms of order 7"(r)/r^"+^, with n > 2. Now, since 7(r) ~ r^'' with 
C < 1, then the second term in ( ^9]) is much smaller than the first one and the omitted terms 
are of higher order. We can therefore keep only the first term and we obtain /(cj) ~ — y | |, 
i.e. the zeroth order approximation for f{ui) is a good one. The validity of this approximation 
of course relies on the fact that C < 1; but, as we shall show below, the result obtained indeed 
coincides with the Imry-Ma argument, i.e. ( = 1/3. If we now substitute this result for f{uj) 
in equation (|i6| ) we can perform the integral (with Qq = 0) and obtain a closed equation for 
q{u). This equation reads 



q{u) 



,-.4 A 2 1/2 \ 1/3 



g2 ^1/2 (-y-) / 4 y/^T^^^Y^ 

— arctanh ' 



4v^ [q]{u) 



-1/3 



(50) 
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where r{u) = — Notice that r{u) is positive at small u since [q]{u) when 

M — s> 0. If r{u) becomes negative, this equation is still valid, provided that we write r^/^(-u) = 
i I r{u) 1^/^ and analytically continue (|50D to imaginary values. This equation relates [q\{u) 
to q{u), with = uq{u) by definition, through a transcendental function and cannot be 

easily solved numerically. It reduces to the equation discussed by Mezard and ParisiB if one 
takes the value of the dipolar interaction to zero (i.e. = 0) and to the equation discussed 
by Hazareesing and Mezardi if one takes the line tension to zero (i.e. a = 0). In these two 
cases, this equation reduces to trivial algebraic equations, which can be easily solved. In 
the first case, one obtains C = 1) which coincides with the exact result of Zhangi and in the 
second case one obtains C = 1/3? as pointed out above. In our case, we can nevertheless 
obtain information regarding the behaviour of q{u) at small u directly from (|50D . Performing 
some algebraic manipulations with (^) and using the identity arctanh(x) = log (j^), one 
can show that 



limig](t.)g-^(«) = -3V^, (51) 



which shows that, since [q\{u) when m — >■ 0, then one must have q{u) = Au^^'^ + 0(m^/^) 
for small u where 0(u^/^) indicates terms of order higher than 1/2 and A — ' — 



12{fcflT)3^ 

Likewise, one must have [q]{u) = ^Au^^"^ + O^u^^"^), where 0(m^/^) indicates terms of order 
higher than 3/2. This coincides with the result obtained by Hazareesing and Mezardi, 
although their solution is valid in a finite interval around 0. In order to obtain a differential 
equation valid Vu G [0, 1], we differentiate equation (|50|) with respect to u. We get 
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q'{u) = 



(52) 



or 



r{u) ~ ^ 2^ ^^^^ 



+ 



12 ^mu) V {ksTr 



g2 ^l/2(-^^ 



4v^ [q]{u) 



arctanh 



4v^ 



u) a 



-4/3 



Now, the second equation can, with the aid of (0), be written in the following form 



[qf{u) = V \q\^[u) + -U q{u) ( [q]\u) - t] q\u) ) 



(53) 



with u = and rj = \^^r^- This equation is also valid when r{u) < without the need 
to perform an analytic continuation, which is an advantage with regard to (0). Now, the 
physical dimensions of v and rj are, according to these definitions, given by [v] = Jm~^ and 
[f]] = J~^m^ (in MKS units). Furthermore, since 5*0 has the dimensions of an energy and 
both x{t) and r have the dimensions of a length, it is also easy to see that the dimensions 
of q{u) are the same as those of u. If we define z{u) = dvq{v), one has q{u) = ^ and 



[q]{u) = u -^ — z{u). Substituting these equations in (|53D , one obtains a differential algebraic 



equation for z{u). This equation reads 



/ dz 



\ -^] + \ u — - z 

\du J \ du 



dz 



dz 

u — -2{z + u) 
du 



0. 



(54) 



Since z{u) and v have the same dimensions, it is useful to write z{u) in terms of a dimension- 
less function z{u) = v x"^ /(x), where x = |^t7¥ is a scahng variable with a = rju = gg 
being a dimensionless constantil. In terms of l{x), the condition ([5l| ) becomes Z(0) = 1/2^/^, 
since the power corresponds to a u^^^ behaviour at small u. Substituting z(u) in terms 



of l{x) in (|54D we obtain the following differential algebraic equation for l{x) 

{xl'{x)+2l{x)Y+ (^xl'{x) + h{x)^ (^x^l'{x)-^x^l{x) -2 1= 



(55) 
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which does not contain any dimensional parameters. The advantage of this equation with 
regard to (^) or ( |54D is its scahng form, which means that once we have obtained the scahng 



function l{x), one simply has to substitute x by its expression in terms of u to obtain z{u) 
for all I' and a. Furthermore, since we have to solve for u G [0,1], this means that 
we need to solve (^) for x G [0, ^^1/4 ]- When a is large0, this interval becomes a small 
interval around and we can linearize (|55D , by writing l{x) = + m(x). Substituting this 
equation above and keeping only terms linear in m{x) and m'(x), we obtain the differential 
equation 

+ — — m (x) + I — -= — —— m{x) = -= , (56j 



8V2 512 J \4V2 256 ; 1024^2 



with solution 



85x2 - 128x^ logx - 3^2: 



X 



4 



"^^^^ 8(3272 + 3x2)2 ' ^^^^ 

where B is an arbitrary constant. Since the solution q'{u) = can be valid in a subinterval 
of [0, 1], the general solution of (|^) is 

z'{u) if u < Uc 

q{u) = { (58) 

z'(Ur) if U > Ur 



where z{u) = v x l{x) with l(x) being given by the full solution of (|55D in the general 
case or simply by the linear approximation l{x) = + ^^^^ ~sli2%^+3x^)^^^'^ ^^^^ 
large a. The values of the arbitrary integration constant B appearing in l{x) and of Uc can 
be determined by substituting the expression for q{u) obtained from z{u) in the equations 
(^) and (|50|). We have solved equation ( p5D numerically using the MANPAK algorithm^, 
available at netlib.org. A comparison between a poljTiomial fit of the numerical solution and 
the linearized solution (|57D is shown in figures 1 and 2, which shows the correctness of this 

23 



solution at small x. However, even in the case of large a, where the linear approximation 
can be used, we still need to solve the system composed by the transcendental equations 
(^) and (|50|), and one would have to resort to a numerical approach. Nevertheless, since 
the behaviour of 7(r) at large length scales is solely determined by the behaviour of [q\{u) at 
small M, it is still possible to show, following Hazareesing and Mezardi, that the roughening 
exponent C, = 1/3. In order to do that, one has to notice that [q]{u) = l/2vx^l{x) + 
3/4z/x'^/'(x) for u < Uc, [q]iu) = [q]{uc) for u > Uc- Substituting this result for [q]{u) in 
equation (^4]), we obtain for small tu, using /(O) = 1/2^/2 g^n^ //(^o) = 0, which follow from 
the linear approximation, the result 

g{uj) = C[ auo^ + (eV2) | uo \ Y'''^ + . . . (59) 

where the dots indicate terms which diverge less strongly at small u) and where C is a 
numerical constant. From this result and from the definition of 7(t), one immediately 
concludes that C = 1/3, in agreement with the Imry-Ma argument of Robbins and Joanny. 
The determination of the crossover behaviour of 7(t), at length scales comparable with 
^ as well as the determination of ^ itself (the Imry-Ma argument just gives the order of 
magnitude) requires, as stated above, the solution of the system of transcendental equations 
(H) and (11). 

VI. CONCLUSIONS 

We have obtained an expression for the elastic energy of a Neel domain-wall in a thin fer- 
romagnetic film in the presence of dipolar interactions and a quenched random field, beyond 
the quadratic approximation for the dipolar energy. Using the replica trick and a varia- 
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tional ansatz, we have obtained a set of self-consistent equations for the Green's functions of 
the displacement field of the domain wall. These equations were solved analytically in the 
case of two-dimensional dipolar interactions by making a quadratic approximation for the 
dipolar energy, which was justified on the basis of the different behavior of the domain- wall 
at different length scales. The problem is then analogous to the one of a domain wall in 
a ferroelastic as studied by Kolomeiskii et al. and to the one of a liquid-gas contact line 
in a disordered solid, as studied by Hazareesing and Mezard. We therefore obtain a value 
( = 1/3 for the roughening exponent of the domain- wall. 

From these calculations, we have obtained some important results. Firstly, we were able 
to represent the dipolar interaction in the domain- wall Hamiltonian as the interaction of the 
quantum particle with a annealed gauge field. Also, in 1-1-1 dimensions, one can integrate 
out such a field, leaving us with polaron like quantum Hamiltonian. Secondly, we have shown 
that a generalized replica symmetry breaking ansatz allows for the treatment of this prob- 
lem. Note that this generalization leads to relatively simple equations in our case, since the 
non-diagonal part of the self-energy in replica space is still wave-vector independent which 
is due to the short-range correlated nature of the random field. The treatment of long-range 
correlated disorder is, on the other hand, a much more complicated problem!!. Finally, we 
were able to justify that taking the quadratic approximation for two-dimensional dipolar 
interactions is sufficient to obtain the correct value of the roughening exponent within this 
approximation. 
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APPENDIX: COUPLING CONSTANTS AND THEIR PHYSICAL UNITS 



Here we collect the definitions of the different couphng constants used in this paper, 
together with their dimensions in the MKS system, denoted by [ ]. One uses the units 
Joule (J=Kgm^s^^), Ampere (A) and meter (m). The 'fundamental' physical constants are 
the magnetic permitivity of the vacuum fiQ, the Bohr magneton fiB, the Lande factor for 
the system ql, the numerical value of the spin S, the strength of the exchange interaction 
between nearest neighbours J (given in J), the amplitude of the random field correlations 
A (given in A^m~^) and the lattice constant a. Furthermore, in the case of a SC/FM/SC 
heterostructure, one also needs the London penetration depth of the superconductor Xl 
and its relative permitivity ^r- From the units of these quantities in the MKSi, one then 
determines the dimension of the coupling constants appearing below. 

a = ^-^^ [a] = Jm-\ (Al) 
a 

(X 



2 fir Xl Oj'^ 



eL= ^":."?' [eL]=Jm-, (A3) 



eL= ^°^\f/^' [eL] = Jm-\ (A4) 
" - ''^'i ^v] = Jm-3, (A5) 



16 a 



[rj] = J-^m^, (A6) 



' 6 /i4 A2 a 

where e^^ is the value of the coupling constant for two-dimensional dipolar interactions (the 
case treated above) and Cg^ is the value for three dimensional dipolar interactions. 
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Figure Captions 

Figure 1. Linearized (continous line) and numerical (dashed line) solutions of equation ^ 



for the same initial (and arbitrary) condition Z(l) = 0.18095. In the case of the numerical 
solution, the MANPAK algorithm was iterated down to x = 0. It is seen that the two 
solutions converge to one another close to a; = 0. 



Figure 2. Left hand side of equation |5^ evaluated for the linearized (continous line) and 
numerical (dashed line) solutions. It is seen that the linearized solution gives better results 
close to X = but that it deviates significantly for larger values. The numerical solution is 
seen to oscilate around zero. Note that the scale of the plot is 10~^. 
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